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We study the magnetization plateau at a third of the saturation magnetization 
of the 5=1/2 trimerized XX Z spin chain at T — 0. The appearance of the 
plateau depends on the values of the XX Z anisotropy and the magnitude of the 
trimerization. This plateauful-plateauless transition is a quantum phase tran- 
sition of the Berezinskii-Kosterlitz-Thouless type, which is difficult to precisely 
detect from the numerical data. To determine the phase boundary line of this 
transition precisely we use the level crossing of low-lying excitations obtained 
from the numerical diagonalization. We also discuss the 5=1/2 ferromagnetic- 
ferromagnetic-antiferromagnetic chain. 
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We study the magnetization plateau at M = M s /3 (M s is the saturation mag- 
netization) of the 5 = 1/2 trimerized XX Z spin chain described by [1] 
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where — l/2<t< lis the trimerization parameter. The bosonized expression 
of the Hamiltonian (4) has the sine-Gordon form 



where v s is the spin wave velocity of the system, II is the momentum density 
conjugate to <fi, [4>{x), n(x')] = iS(x — x'), and the coefficients v s , K, and are 
smooth functions of J , t and A The field <j) is related to the fast varying (in 
space) part of the spin density S z (x) in the continuum picture as 



which makes it clear the physical meaning of (f>. 

As is well known, the excitation spectrum of the sine-Gordon model is ei- 
ther massive (plateauful) or massless (plateauless) depending on the values of 
parameters. These two states are distinguished by the renormalized value of 
K: the former is realized when K < 4 and the latter K > 4. The transition 
between these two states is of the Brerezinskii-Kosterlitz-Thouless type, which 
was first pointed out by one of the present authors (K.O) [2]. To determine 
the K = 4 BKT transition point from the numerical diagonalization data for 
finite systems, we use the method developed by Nomura and Kitazawa [3]. They 
discussed the K = 4 BKT transition at M = to find that the crossing be- 
tween the M = excitation with the twisted boundary condition (TBC) and 
the M = 2 excitation with the periodic boundary condition (PBC). Because we 
discuss the finite magnetization case now, it is necessary to use the Legendre 
transformation E — > E— HM. Thus we can conclude that the crossing between 
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SLt(x) = - I cos (2k F x - | + V2V) + ^ 
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where E% BC = (l/2){£j BC (M s /3 + 2) + £^ BC (M s /3 - 2)}. Figure 1 shows 
the crossing between these two excitations, from which we see A C (L = 18) = 
—0.8389. By sweeping parameters, we obtain the phase diagram on the A — t 
plane (Fig.l). 
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Figure 1 : Ai? TBC (opne circles) 
and AE PBC (closed circles) for 
L = 18 spins as functions of 
anisotropy parameter A when 
t = —0.25. From the crossing 
point we obtain A C (L = 18) = 
-0.8389. 



Figure 2: Phase diagram on 
the A — t plane. Closed circles 
are the BKT transition points 
determined from the numerical 
data as explained in the text. 
Thick dotted line denotes the 
Gaussian line. The multicrit- 
ical point M is at (A, t) = 
(-0.729,0). 



We also investigate the ferromagnetic-ferromagnetic-antiferromagnetic model 
[2,4] 



H F FA = ^{-J F [/l3j-2,3j-l(l) + /i3j-l,3j(l)] 
3 

+ ^A./i3j,3j+l(l)} 



(9) 



where Jp and J a are the ferromagnetic and antifcrromagnctic interaction con- 
stants, respectively. We set A = 1 so that -Hffa is reduced to the isotropic 
5 = 3/2 antiferromagnetic chain when 7 = Jf/Ja — ► 00. Following Hida's nu- 
merical calculation, the M s /3 plateau clearly exists when 7 = Jf/Ja is small. 
On the other hand, it is believed that there exists no M s /3 plateau in the 
5 = 3/2 chain. Thus the plateauful-plateauless transition takes place at finite 
7. Since the Hamiltonian Hffa is transformed into the generalized version of the 
Hamiltonian (1) by the spin rotation [2], the present method is also applicable 
to Hffa, resulting in A c = 15.4 [5]. 

In conclusion, we have obtained the phase diagram for Hamiltonian (1), 
developing a new method to detect the plateauful-plateualess transition from 
the numerical diagonalization data. Our new method is applicable to various 
systems. 
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